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ABSTRACT 



We present a measurement of the 6-quark inclusive fragmentation function in Z decays using a novel 
kinematic I?-hadron energy reconstruction technique. The measurement was performed using 350,000 
■ hadronic Z° events recorded in the SLD experiment at SLAC between 1997 and 1998. The small and 
stable SLC beam spot and the CCD-based vertex detector were used to reconstruct I?-decay vertices 
with high efficiency and purity, and to provide precise measurements of the kinematic quantities used in 
this technique. We measured the B energy with good efficiency and resolution over the full kinematic 
range. We compared the scaled i?-hadron energy distribution with models of 6-quark fragmentation 
o 1 and with several ad hoc functional forms. A number of models and functions are excluded by the data. 
The average scaled energy of weakly-decaying B hadrons was measured to be < Xb > = 0.709 ± 0.003 
(stat) ± 0.003 (syst) ± 0.002 (model). 
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1 Introduction 



The production of heavy hadrons (H) in e + e~ annihilation provides a laboratory for the study of 
heavy-quark (Q) jet fragmentation. This is commonly characterised in terms of the observable xjj = 
IE n j y/s, where Ejj is the energy of a B or D hadron containing a b or c quark, respectively, and \/s 
is the cm. energy. In contrast to light-quark jet fragmentation one expects JTJ the distribution of xh, 
D(xh), to peak at an x^-value significantly above 0. Since the hadronisation process is intrinsically 
non-perturbative D(xh) cannot be calculated directly using perturbative Quantum Chromodynam- 
ics (QCD). However, the distribution of the closely- related variable xq = 2Eq/^/s can be calculated 
perturbatively ||, [3|, |J ]5| and related, via model-dependent assumptions, to the observable quan- 
tity D(xh)', a number of such models of heavy-quark fragmentation have been proposed || [7|, ||, |9|. 
Measurements of D(xh) thus serve to constrain both perturbative QCD and the model predictions. 
Furthermore, the measurement of D(xh) at different cm. energies can be used to test QCD evolution, 
and comparison of D(xb) with D(xd) can be used to test heavy-quark symmetry ||10|| . Finally, the 
uncertainty on the forms of D(xd) and D(xb) must be taken into account in studies of the produc- 
tion and decay of heavy quarks, see eg. [p]] ; more accurate measurements of these forms will allow 
increased precision in tests of the electroweak heavy-quark sector. 

We have measured the inclusive weakly-decaying 5-hadron scaled energy distribution D(xb) in Z° 
decays. Earlier studies 1^] used the momentum spectrum of the lepton from semi-leptonic B decays 
to constrain the mean value < xb > and found it to be approximately 0.70; this is in agreement with 
the results of similar studies at \fs = 29 and 35 GeV fl3| . In more recent analyses |14|, [15], [L6|] D(xb) 
has been measured by reconstructing B hadrons via their B — > DIX decay mode. In this case the 
reconstruction efficiency is intrinsically low due to the small branching ratio for B hadrons to decay 
into the high-momentum leptons used in the tag. Also, the reconstruction of the 5-hadron energy 
using calorimeter information usually has poor resolution for low B energy, resulting in poor sensitivity 
to the shape of the distribution in this region. 

We present the results of a new method for reconstructing I?-hadron decays and the B energy 
inclusively, using only charged tracks, in the SLD experiment at SLAC. We used the upgraded charge- 
coupled device (CCD) vertex detector, installed in 1996, to reconstruct S-decay vertices with high 
efficiency and purity. Combined with the micron-size SLC interaction point (IP), precise vertexing 
allowed us to reconstruct accurately the B flight direction and hence the transverse momentum of 
tracks associated with the vertex with respect to this direction. Using the transverse momentum and 
the total invariant mass of the associated tracks, an upper limit on the mass of the missing particles 
was found for each reconstructed S-decay vertex, and was used to solve for the longitudinal momentum 
of the missing particles, and hence for the energy of the B hadron. In order to improve the B sample 
purity and the reconstructed 5-hadron energy resolution, B vertices with low missing mass were 
selected. The method is described in Section 3. In Section 4 we compare our reconstructed D(xb) 
with the predictions of heavy-quark fragmentation models. We also test several functional forms for 
this distribution. In Section 5 we describe the unfolding procedure used to derive our estimate of the 
true underlying D(xb)- In Section 6 we discuss the systematic errors. In Section 7 we summarize the 
results. Our measurement based on a data sample one-third the size of that used here is reported in 
Refs. 1X71, 
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2 Apparatus and Hadronic Event Selection 



This analysis is based on roughly 350,000 hadronic events produced in e + e~ annihilations at a mean 
center-of-mass energy of yfs = 91.28 GeV at the SLAC Linear Collider (SLC), and recorded in the SLC 
Large Detector (SLD) in 1997 and 1998. A general description of the SLD can be found elsewhere |19[ 



The trigger and initial selection criteria for hadronic Z° decays are described in Ref. ||20|| . This analysis 
used charged tracks measured in the Central Drift Chamber (CDC) |^TJ and in the upgraded Vertex 
Detector (VXD3) f22"| . Momentum measurement was provided by a uniform axial magnetic field of 
0.6T. The CDC and VXD3 give a momentum resolution of cr p± /p± = 0.01 © 0.0026p^, where p± is 
the track momentum transverse to the beam axis in GeV/c. In the plane normal to the beamline 
the centroid of the micron-sized SLC interaction point (IP) was reconstructed from tracks in sets of 
approximately thirty sequential hadronic Z° decays with a precision of a r f p ~ 4 ± 2 /im. The IP 
position along the beam axis was determined event by event using charged tracks with a resolution of 
cxfp ~ 20 /mi. Including the uncertainty on the IP position, the resolution on the charged-track impact 
parameter (d) projected in the plane perpendicular to the beamline was = 80)33 / (p sm 3 ^ 2 9) /im, 
and the resolution in the plane containing the beam axis was o z d = 10©33/(psin 3 / 2 9) /im, where 9 is 
the track polar angle with respect to the beamline. The event thrust axis p3| was calculated using 



energy clusters measured in the Liquid Argon Calorimeter [f2~4]j . 

A set of cuts was applied to the data to select well-measured tracks and events well contained 
within the detector acceptance. Charged tracks were required to have a distance of closest approach 
transverse to the beam axis within 5 cm, and within 10 cm along the axis from the measured IP, as 
well as | cos#| < 0.80, and p± > 0.15 GeV/c. Events were required to have a minimum of seven such 
tracks, a thrust-axis polar angle w.r.t. the beamline, 9t, within | cos9t\ < 0.71, and a charged visible 
energy E V i S of at least 20 GeV, which was calculated from the selected tracks, which were assigned 
the charged pion mass. The efficiency for selecting a well-contained Z° —>■ qq(g) event was estimated 
to be above 96% independent of quark flavor. The selected sample comprised 218,953 events, with an 
estimated 0.10 ± 0.05% background contribution dominated by Z° — > t + t~ events. 

For the purpose of estimating the efficiency and purity of the 5-hadron selection procedure we 
made use of a detailed Monte Carlo (MC) simulation of the detector. The JETSET 7.4 event 
generator was used, with parameter values tuned to hadronic e + e~ annihilation data p6| , combined 
with a simulation of B hadron decays tuned |^7| to T(4S) data and a simulation of the SLD based on 



GEANT 3.21 p8[ . Inclusive distributions of single-particle and event-topology observables in hadronic 
events were found to be well described by the simulation [pOj . Uncertainties in the simulation were 
taken into account in the systematic errors (Section |6|). 



3 5-Hadron Selection and Energy Measurement 
3.1 B-Hadron Selection 

The B sample for this analysis was selected using a 'topological vertexing' technique based on the 
detection and measurement of charged tracks, which is described in detail in Ref. p9] . Each hadronic 
event was divided into two hemispheres by a plane perpendicular to the thrust axis. In each hemisphere 
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the vertexing algorithm was applied to the set of 'quality' tracks having (i) at least 23 hits in the CDC 
and 2 hits in VXD3; (ii) a combined CDC and VXD3 track fit quality of xV^W <8; (iii) a momentum 
(jp) in the range 0.25< p <55 GeV/c, (iv) an impact parameter projection in the r — plane of less 
than 0.3 cm, and a projection along the z axis of less than 1.5 cm; (v) an r — impact parameter 
error no larger than 250 /xm. 

Vertices consistent with photon conversions or K° and A decays were discarded. In hemispheres 
containing at least one found vertex the vertex furthest from the IP was retained as the 'seed' vertex. 
Those events were retained which contained a seed vertex separated from the IP by between 0.1 cm 
and 2.3 cm. The lower bound reduces contamination from non-P-decay tracks and backgrounds from 
light-flavor events, and the upper bound reduces the background from particle interactions with the 
beam pipe. A sample of 76,421 event hemispheres was selected. 

In each hemisphere, a vertex axis was defined as the straight line joining the IP to the vertex, which 
was located at a distance D from the IP. For each quality track not directly associated with the vertex, 
the distance of closest approach to the vertex axis, T, and the distance from the IP along the vertex 
axis to the point of closest approach, L, were calculated. Tracks satisfying T < 1 mm and L/D > 0.3 
were added to the vertex. These T and L cuts were chosen to minimize false track associations to 
the seed vertex, since typically the addition of a false track has a much greater kinematic effect than 
the omission of a genuine P-decay track, and hence has more effect on the reconstructed P-hadron 
energy. Our Monte Carlo studies show that, on average, this procedure attaches 0.85 tracks to each 
seed vertex, 91.9% of the tracks from tagged true P decays are associated with the resulting vertices, 
and 98.0% of the vertex tracks are from true P decays. 

The large masses of the B hadrons relative to light-flavor hadrons make it possible to distinguish 
P-hadron decay vertices from those vertices found in events of primary light flavor using the vertex 
invariant mass, M. However, due to the effect of those particles missed from being associated with the 
vertex, which are mainly neutrals, M cannot be fully determined. In the rest frame of the decaying B 
hadron, M can be written 



where M ch and M are the total invariant masses of the set of vertex-associated tracks and the set of 
missed particles, respectively. P t is the momentum sum, transverse to the B flight direction, of the 
vertex-associated tracks, which, by momentum conservation, is identical to the transverse momentum 
sum of the missed particles. P c hi and Pqz are t ne respective momentum sums along the B flight 
direction. In the B rest frame, P c u = P z- Using the set of vertex-associated charged tracks, we 

— * 

calculated the total momentum vector P ch and its component transverse to the flight direction P t , and 
the total energy E ch and invariant mass M ch , assuming the charged-pion mass for each track. The 
lower bound for the mass of the decaying hadron, the 'P r corrected vertex mass', 



was used as the variable for selecting B hadrons. Our simulations show that the majority of non-P 
vertices have Mp t less than 2.0 GeV/c 2 . However, occasionally the measured P t may fluctuate to a 
much larger value than the true P t , causing some charm-decay vertices to have Mp t larger than 2.0 




(1) 




(2) 
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GeV/c 2 . To reduce this contamination, we calculated the 'minimum Pt by allowing the IP and the 
vertex to float to any pair of locations within the respective one-sigma error ellipsoids. We substituted 



the minimum P t in Equation (fj) and used this modified Mp t as our variable for selecting B hadrons [30 

Figure [l] shows the distribution of Mp t for the selected sample of hemispheres containing a vertex, 
and the corresponding simulated distribution. P-hadron candidates were selected by requiring Mp t 
> 2.0 GeV/c 2 . We further required Mp t < 2 x M c h to reduce the contamination from fake vertices 
in light-quark events |30|| . A total of 42,093 hemispheres were selected, with an estimated efficiency 



for selecting a true 5-hemisphere of 43.7%, and a sample purity of 98.2%. The contributions from 
light-flavors in the sample were 0.34% for primary u, d and s hemispheres and 1.47% for c hemispheres. 

3.2 B-Hadron Energy Measurement 

The energy of each B hadron, Ep, can be expressed as the sum of the reconstructed- vertex energy, 
E c h, and the energy of those true 5-decay particles that were missed from the vertex, Eq. Eq can be 
written 

El = M 2 + P 2 + P 2 (3) 

The two unknowns, M and Pqi, must be found in order to obtain E . One kinematic constraint can 
be obtained by imposing the P-hadron mass, Mp, on the vertex. From Equation ([I]) we derive the 
following inequality, 



>Ml h + Pf + ^M 2 + Pi < M B , (4) 

where equality holds in the limit that both Pqi and P c u vanish in the P-hadron rest frame. Equation (|J) 
effectively sets an upper bound on M , M 0max : 



M% max = M%- 2M BV /Mi + Ff + M 2 ch . (5) 
The lower bound is zero. Hence 

< Ml < M% max , (6) 

f 2 

Omax 



and we expect to obtain a good estimate of M , and therefore of the P-hadron energy, when MQ max is 



small. 

We used our simulation to study this issue. We find that the true value of M tends to cluster near 
its maximum value Mo max . Figure |2| shows the relative deviation of Mo max from the true Mq for all 
B hadrons, assuming Mp = 5.28 GeV/c 2 in Eq. Although approximately 20% of the B hadrons 
are B Q S and A&, which have larger masses than the B° and B^, the values of M 0max obtained using 
Afg=5.28 GeV/c 2 are typically within about 10% of M . The distribution of the reconstructed Ml max 
for vertices in the selected-hemisphere sample is shown in Figure the negative tail is an effect of 
detector resolution. The simulation is in good agreement with the data, and implies that the non-Z? 
background is concentrated at high Ml max ; this is because most of the light-flavor vertices have small 
Mp t and therefore, due to the strong negative correlation between Mp t and M Qmax , large M 0max . 

Because, for true B decays, M peaks near M Qmax , we set Mq = Ml max if Ml max >0, and Mq = 
if Ml max <0. We then calculated P f. 

p _ jg ~ (jg> + Pt 2 ) ~ (Mq 2 + PI) p m 
01 ~ 2(M 2 , + P 2 ) chh [) 
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and hence E (Equation (||)). We divided the reconstructed 5-hadron energy, E r B ec = E + E c h, by the 
beam energy, E beam , to obtain the reconstructed scaled 5-hadron energy, x r B c . 

The resolution of x r B c depends on both M^ max and the true xb, x% ue . Using our simulation we 
found that vertices that have Ml max < — 1.0(GeV/c 2 ) 2 are often poorly reconstructed; we rejected 
them from further analysis. Vertices with small values of |Mg ma J are typically reconstructed with 
better resolution and an upper cut on MQ max was hence applied. For an x^-independent M^ max cut 
we found that the efficiency for selecting B hadrons is roughly linear in i^- in order to obtain an 
approximately x^-independent selection efficiency we required: 

Km* < {1 1 + 0.007(£ beam - E r B ec ) + A.0exp[-(E r B ec - 5.5)/3.5]} 2 , (8) 

where the two ad hoc terms that depend on E B ec increase the efficiency at lower S-hadron energy. 

In addition, in order to reduce the light-flavor background, each vertex was required to contain 
at least 3 quality tracks with a normalized impact parameter greater than 2. This cut reduces the 
dependence of the reconstructed I?-hadron energy distribution on the light-flavor simulation in the 
low-energy region. 

A total of 4,164 hemispheres contained vertices that satisfied these selection cuts. Figure [| shows 
the distribution of Mg maa: ; the simulation and data are in good agreement. We calculated that the 
efficiency for selecting B hadrons is 4.17% and the I?-hadron purity is 99.0%, with a uds (charm) 
background of 0.4% (0.6%). The efficiency as a function of the true xb value, x l B ue , is shown in 
Figure |5|. The dependence is weak except for the lowest xb region; the efficiency is substantial, even 
just above the kinematic threshold. 

We examined the energy resolution of this technique using simulated events. The distribution of the 
normalized difference between the true and reconstructed scaled 5-hadron energies, (x r B c — x% ue )/x% ue , 
was fitted with the sum of two Gaussians. A feature of the analysis is that the distribution is symmetric 
and the fitted means are consistent with zero. The fit yields a core width (the width of the narrower 
Gaussian) of 9.6% and a tail width (the width of the wider Gaussian) of 21.2%, with the narrower 
Gaussian representing a population fraction of 83.6%. Figure || shows the core and tail widths as a 
function of x l B ue , where, in order to compare the widths from different xb bins, the ratio between the 
core and tail populations was fixed to that obtained above. The x^-dependence of the resolution is 
weak. The resolution is good even at low B energy, which is an advantage of this energy reconstruction 
technique. 

Figure [7| shows the distribution of the reconstructed scaled 5-hadron energy; the simulated dis- 
tribution is also shown. The small non-5 background, the high B selection efficiency over the full 
kinematic coverage, and the good energy resolution combine to give a much improved sensitivity of 
the data to the underlying true shape of the B energy distribution (see next section). The distribution 
of the non-i? background was subtracted bin-by-bin to yield D rec (x r B c ), which is shown in Fig. |8]. 

The JETSET event generator used in our simulation is based on a perturbative QCD 'parton 
shower' for production of quarks and gluons, together with the phenomenological Peterson function || 
(Table |TJ)0 to account for the fragmentation of b and c quarks into B and D hadrons, respectively, 
within the iterative Lund string hadronisation mechanism p3| . It is apparent that this simulation 

*We used a value of the Peterson function parameter ej, = 0.006 Q. 
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does not reproduce the data (Figure |7|); the y 2 for the comparison is 70.3 for 16 binJ. 



4 The Shape of the 5-Hadron Energy Distribution 
4.1 Tests of 6-Quark Fragmentation Models f(z,/3) 

We tested models of 6-quark fragmentation. Since the resulting fragmentation functions are usually 
functions of an experimentally-inaccessible variable z, eg. z — (E + p\\)h/(E + P\\)q or z = P\\ h /p\\q, 
where p\\ represents the hadron momentum along the primary heavy-quark momentum vector, it is 
necessary to use a Monte Carlo generator to produce events according to a given input fragmentation 
function f(z,/3), where (3 represents the set of model arbitrary parameters. 

We considered the phenomenological models of the Lund group ||, Bowler M, Peterson et al. || 
and Kartvelishvili et al. @. We also considered the perturbative QCD calculations of Braaten et 
al. (BCFY) H|, and of Collins and Spiller (CS) @. Table contains a list of the models. We 
implemented in turn each fragmentation model in JETSET and generated events without detector 
simulation. In addition, we tested the UCLA fragmentation model 3S| with default parameter settings, 



as there is no explicit parameter for controlling the 5-hadron energy. We also tested the HERWIG |33 
event generator, and used both possible settings of the parameter switch cldir. cldir=l forces the 
heavy hadron to continue in the heavy-quark direction in the hadronisation-cluster decay rest frame, 
and thereby hardens the fragmentation function. cldir=0 suppresses this feature and yields a softer 
fragmentation function. 

In order to make a consistent comparison of each model with the data we adopted the following 
procedure. For each model starting values of the arbitrary parameters, /3, were assigned and the cor- 
responding fragmentation function f(z,(3) was used to produce the scaled weakly-decaying S-hadron 
energy distribution, -D^odei ( x s Me , 0) before simulation of the detector. The corresponding reconstructed 
distribution, D^ del (x r ^ c , (3), was derived from the reconstructed distribution generated with our de- 
fault model, D^ ault (x r £ c ) (Fig. |7|), by weighting events at the generator level with the weight factor 
Dmodei(. x B Ue i P) / Dd™f e a uit(, xt B Ue ) ■ The resulting reconstructed distribution was then compared with the 
data distribution, and the \ 2 value, defined as 

data _ JVfMC \ 2 

- (9) 




was calculated, where N is the number of bins used in the comparison, Nf ata is the number of entries in 
bin % in the data distribution, and Nf' IC is the number of entries in bin i in the simulated distribution!. 
(7j is the statistical error on the deviation of the observed number of entries for the data from the 
expected number of entries in bin i, which can be expressed as 

(fNT 5 ) ^+(r V / A^) 2 , (10) 



a 2 



^We excluded from the comparison several bins that contained very few events; see Section 4.1. 
*r is the factor by which the total number of entries in the simulated distribution was scaled to the number of entries 
in the data distribution; r ~ 1/12. 
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Model 


f(z,P) 


Reference 


BCFY 


[l-(lV^]«[ 3 + ^i(-^/i(r)] 


1 


Bowler 




Q 




v 2 1 1 - z 1 - z 1 _ ^ 




Kartvelishvili 


z Qb (l -z) 




Lund 


4(1 — -2) a exp(— hm\jz) 


i 


Peterson 




§ 



Table 1: 6-quark fragmentation models used in comparison with the data. For the BCFY model, 
f^r) = 3(3 - 4r), f 2 (r) = 12 - 23r + 26r 2 , / 3 (r) = (1 - r)(9 - llr + 12r 2 ), and 
/ 4 (V) = 3(l-r) 2 (l-r + r 2 ). 



where \^yrN^ Ic j is the expected statistical variance on the observed number of entries in bin z, 

assuming the model being tested is correct, and (r^N^ } is the statistical variance on the expected 

number of entries in bin i. Since the x 2_ test is not statistically effective for bins with a very small 
number of entries, the third, the fourth, and the last three bins in Figure [7| were excluded from the 
comparison. 

For each model we varied the values of the parameters j3 and repeated the above procedure. The 
minimum y 2 was found by scanning through the input parameter space, yielding a set of parameters 
which give an optimal description of the reconstructed data by the fragmentation model in question. 
The resulting distributions are shown in Figure |8|. Table ^ lists the results of the comparisons. 

We conclude that with our resolution and our current data sample, we are able to distinguish among 
these fragmentation models. Within the context of the JETSET fragmentation scheme, the Lund and 
Bowler models are consistent with the data with x 2 probabilities of 31% and 35%, respectively, the 
Kartvelishvili model is consistent with the data at the 1% level, while the Peterson, BCFY and CS 
models are found to be inconsistent with the data. The UCLA model is consistent with the data at 
a level of 6% \ 2 probability. The HERWIG model with cldir=0 is confirmed to be much too soft; 
using cldir=l results in a harder distribution and a substantial improvement, but it is still too soft 
relative to the data. 
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Model 


X 2 /dof 


Parameters 


(x B ) 


JETSET + BCFY 


105/16 


r = 0.085 


0.694 


JETSET + Bowler* 


17/15 


a = 1.4, b = 1.2, (r 6 = 1) 


0.709 


JETSET + CS 


142/16 


e b = 0.003 


0.691 


JETSET + Kartvelishvili* et al. 


32/16 


a b = 10.0 


0.708 


JETSET + Lund* 


17/15 


a = 1.4,6 = 0.4 


0.712 


JETSET + Peterson et al. 


70/16 


e b = 0.0055 


0.700 


HERWIG cldir=0 


1015/17 




0.632 


HERWIG cldir=l 


149/17 




0.676 


UCLA* 


27/17 




0.718 



Table 2: Results of fragmentation model tests. The minimum x 2 , number of degrees of freedom, 
corresponding parameter values, and the mean value of the corresponding 5-energy distribution are 
listed. A * indicates those models used below to unfold the data. 



4.2 Tests of Functional Forms f(xB,X) 

We considered the more general question of what functional forms, f(xs, A), can be used as estimates 
of the true scaled 5-energy distribution. We considered the functional forms of the BCFY, CS, 
Kartvelishvili, Lund, and Peterson groups in terms of the variable x 1 ^ . In addition we considered 
ad hoc generalisations of the Peterson function ('F'), an 8th-order polynomial ('P8') and a 'power' 
function. These functions are listed in Table 131 Each function vanishes at xb = and Xb — 1- 



Function 



f(x B ,X) 

;i + 6(i-x B )) (1 c d 



Reference 



F 
P8 

Power 



xb y ^ x B 1 - xb' 

X B {1 - X B ){X B - X B )(l + T,i=lPiX B ) 

x%{i-x B y 







(see text) 
(see text) 



Table 3: Additional ad hoc B energy functional forms used in comparison with the data; x° B = 

Mb/ Ebeam- 

For each functional form, a testing procedure similar to that described in subsection [0] was 
applied. The fitted parameters A and the minimum x 2 values are listed in Table |j. The corresponding 
D model (^b °) are compared with the data in Figure |]. 
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Function 


X 2 /dof 


Parameters 


(xb) 


BCFY 


73/16 


r = 0.248 ±0.007 


0.704±0.003 


CS 


75/16 


e b = 0.0519 ±0.0036 


0.706±0.003 


Kartvelishvili et al. 


138/16 


a b = 3.904 ±0.072 


0.710±0.003 


Lund 


252/15 


a = 1.88 ±0.08 
hm\ = 0.32 ±0.05 


0.715±0.003 


Peterson et al* 


31/16 


e b = 0.0382 ±0.0016 


0.709±0.003 


Fl* 


20/15 


c = 0.884 ±0.014 
d = 0.0181 ±0.0015 


0.707±0.003 


F2* 


31/15 


c = 0.976 ±0.029 
d = 0.039 ±0.002 


0.710±0.003 


P8* 


12/12 


Pl = -9.99 ±0.25 
p 2 = 40.84 ± 0.25 

p 3 = -82.26 ± 0.68 
p 4 = 80.90 ± 0.76 

p 5 = -30.60 ± 0.54 


0.709±0.003 


Power 


133/15 


a = 3.73 ±0.17 
P = 0.84 ±0.07 


0.713±0.003 



Table 4: Results of the x 2 fit of ad hoc functions to the reconstructed scaled _B-hadron energy distribu- 
tion. The minimum x 2 value, the number of degrees of freedom, the corresponding parameter values, 
and the mean value of the corresponding 5-energy distribution are listed. Errors are statistical only. 
A * indicates those functions used below to correct the data. 
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Two sets of optimised parameters were found for the generalised Peterson function F: 'FT, obtained 
by setting the parameter b (Table ||) to infinity, behaves like xb as xb — ► and (1 — xb) 3 as xb —* 1 
and yields a y 2 probability of 18%; 'F2', obtained by setting b to zero, has a x 2 probability of 1.0%. 
A constrained polynomial of at least 8th-order was needed to obtain a \ 2 probability greater than 
0.1%. The Peterson function reproduces the data with a x 2 probability of about 1%. The remaining 
functional forms are found to be inconsistent with the data. The widths of the BCFY and CS functions 
are too large to describe the data; the Kartvelishvili, Lund and 'power' functions vanish too fast as 
Xb — > 0. We conclude that, within our resolution and with our current data sample, we are able to 
distinguish among these ad hoc functional forms. 

5 Correction of the 5-Energy Distribution 

In order to compare our results with those from other experiments and potential future theoretical 
predictions it is necessary to correct D rec (x r B C ) for the effects of detector acceptance, event selection 
and analysis bias, as well as for bin-to-bin migrations caused by the finite resolution of the detector 
and the analysis technique. Due to the known rapid variation of the a priori unknown true 5-energy 
distribution at large xb, any correction procedure will necessarily be model-dependent. We chose a 
method that allows explicit evaluation of this model-dependence and which gives a very good estimate 
of the true energy distribution using all of the above models or functional forms that are consistent 
with the data. 

We applied a 25 x 25 matrix unfolding procedure to D rec (x r B C ) to obtain an estimate of the true 
distribution D true (x t B " ae ), where x% ue refers to the weakly-decaying B hadron: 

D true (x% ue ) = e~ 1 {x% ue ) ■ E(x% ue ,x T £ c ) -D rec {x r £ c ) (11) 

where E is a matrix to correct for bin-to-bin migrations, and e is a vector representing the efficiency 
for selecting true S-hadron decays. E and e were calculated from our MC simulation; the matrix E 
incorporates a convolution of the input fragmentation function with the resolution of the detector. 
E(i,j) is the number of vertices with x^ ne in bin i and x r ^ c in bin j, normalized by the total number 
of vertices with x r ^ c in bin j. 

We evaluated E by using in turn the Monte Carlo simulation weighted according to each input 
generator-level true B energy distribution found to be consistent with the data in Section We 
considered in turn each of the eight consistent distributions, using the optimised parameters listed in 
Tables ^| and ^. The matrix E was then evaluated by examining the population migrations of true B 
hadrons between bins of the input scaled B energy, x^ ne , and the reconstructed scaled B energy, x r £ c . 
Using each D^^^x^ 6 ) , the data distribution D rec (x r B C ) was then unfolded according to Equation fllTD 
to yield D true {x t B Ue ) ) which is shown for each input fragmentation function in Figure It can be seen 
that the shapes of D true {x t B Ue ) differ systematically among the input scaled S-energy distributions. 
These differences were used to assign systematic errors. 
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6 Systematic Errors 



We considered sources of systematic uncertainty that potentially affect our measurement. These may 
be divided into uncertainties in modelling the detector and uncertainties on experimental measurements 
serving as input parameters to the underlying physics modelling. For each source of systematic error, 
the Monte Carlo distribution D^^^x 1 ^) was re-weighted and then the resulting new reconstructed 
distribution, D^^x 7 ^ ) , was compared with the data D rec (x r £ c ) by repeating the fitting and unfolding 
procedures described in Sections 4 and 5. The differences in both the shape and the mean value of the 
x true distribution relative to the default procedure were considered. 

Ad hoc corrections were applied to the simulations of four track-related quantities to account for 
discrepancies w.r.t. the data, namely the tracking efficiency and the distributions of track p±, polar 
angle and the projection of the impact parameter along the z axis. In each case a systematic error 
was assigned (see Table ^) using half the difference between the results obtained with the default and 
corrected simulations. 

A large number of measured quantities relating to the production and decay of charm and bottom 
hadrons are used as input to our simulation. In bb events we considered the uncertainties on: the 
branching fraction for Z° — > bb; the rates of production of B^, B° and B® mesons, and b baryons; 
the lifetimes of B mesons and baryons; and the average B hadron decay charged multiplicity. In 
cc events we considered the uncertainties on: the branching fraction for Z° — > cc, the charmed hadron 
lifetimes, the charged multiplicity of charmed hadron decays, the production of K° from charmed 
hadron decays, and the fraction of charmed hadron decays containing no 7r°s. We also considered the 
rate of production of ss in the jet fragmentation process, and the production of secondary bb and 
cc from gluon splitting. The world-average values and their respective uncertainties [|TT| |31| were used 
in our simulation and are listed in Table [|. Most of these variations affect the normalisation, but have 
very little effect on the shape or the mean value. In no case do we find a variation that changes our 
conclusion about which models and functions are consistent with the data. The systematic errors on 
the mean value are listed in Table |5]. 
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Source 


Variation 


5 (x B ) 


Tracking efficiency correction 
Impact parameter smearing in z 
Track polar angle smearing 
Track l/p± smearing 


-1.5±0.75% 
9.0±4.5/im 
1.0±0.5 mrad 
0.8±0.4 MeV" 1 


0.0007 
0.0006 
0.0002 
0.0013 


Detector total 




0.0016 


B~^~ production fraction 
production fraction 

JD S piUClUCLlUIl ilcxCllOIl 
A 7 nrnnnrtinri TT£i r^i/"!!"! 

li-Q \-Ji- vj LI LIKj Ivjll 11 CIKj KjLyJLi. 

B — > charm multiplicity and species 
B — > K° multiplicity 
B — > A multiplicity 

B decay < n c h > 

D — > K° multiplicity 

D — > no 7T° fraction 

D decay < n ch > 

g^bb 

g->cc 

B° mass 

b, c hadron lifetimes, Rf,, R c 


U.oy±U.ll 

n ^Q-i-n 1 1 
u.oyztu. 1 1 

n nos-i-n nni 9 
n i nq-un ni <a 

U. 1UOXU.U X O 

[18| 
0.658±0.066 
0.124±0.008 
4.955±0.062 

[18| 

PI 

18 

0.00254±0.00050 /evt 

0.0299±0.0039 /evt 
5.2794± 0.0005 GeV/c 2 




^U.UUUl 
^■n nnm 

<U.UUU1 

_i_n nnn^ 
0001 

±0.0006 
±0.0009 
±0.0002 

-0.0004 
+0.0001 

±0.0014 
±0.0006 
±0.0003 
±0.0001 
±0.0003 
<0.0001 
±0.0002 


Physics total 




0.0020 


Monte Carlo statistics 




0.0008 


Total systematic 




0.0027 



Table 5: Uncertainty source, range of variation and size of the resulting systematic error on < xb >■ 
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Other relevant systematic effects such as variation of the event selection cuts and the assumed B- 
hadron mass were also found to be very small. As a cross-check, we varied the M 0max cut (Equation (§)) 
used to select the final B sample and repeated the analysis procedure. In each case, conclusions about 
the shape of the B energy distribution hold. In each bin, all sources of systematic uncertainty were 
added in quadrature to obtain the total systematic error. 

The model-dependence of the unfolding procedure was estimated by considering the envelope of 



the unfolded results shown in Figure [K| Since eight models or functions are consistent with the data, 
in each bin of x^ ue we calculated the average value of these eight unfolded results as well as the r.m.s. 
deviation; the average was taken as our central value and the deviation was assigned as the unfolding 
uncertainty. Figure |TT] shows the final corrected xb distribution D(xb)- The data are listed in Table ^ 
Since two of the eight functions (the Kartvelishvili model and the Peterson functional form) are only 
in marginal agreement with the data, and the 8th-order polynomial has an unphysical behavior near 
Xb = 1, this r.m.s. may be considered to be a rather reasonable envelope within which the true xb 
distribution is most likely to vary. The model dependence of this analysis is significantly smaller than 
that of previous direct S-energy measurements, indicating the enhanced sensitivity of our data to the 
underlying true energy distribution. 
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x B range 


D(x B ) 


stat. 


systematic 


unioldmg 


total 


0.00 


< 


x B 


< 


0.04 


0.000 


0.000 


0.000 


0.000 


U.UUU 


0.04 


< 


x b 


< 


0.08 


0.000 


0.000 


0.000 


0.000 


U.000 


0.08 


< 


x b 


< 


0.12 


0.000 


0.000 


0.000 


0.000 


0.000 


0.12 


< 


x b 


< 


0.16 


0.110 


0.029 


0.004 


0.014 


0.034 


0.16 


< 


x b 


< 


0.20 


0.188 


0.035 


0.005 


0.025 


0.043 


0.20 


< 


x b 


< 


0.24 


0.204 


0.032 


0.006 


0.013 


0.036 


0.24 


< 


x b 


< 


0.28 


0.213 


0.027 


0.008 


0.010 


0.030 


0.28 


< 


x b 


< 


0.32 


0.268 


0.031 


0.009 


0.015 


0.036 


0.32 


< 


x b 


< 


0.36 


0.340 


0.036 


0.011 


0.011 


0.039 


0.36 


< 


x b 


< 


0.40 


0.398 


0.037 


0.012 


0.010 


0.041 


0.40 


< 


x b 


< 


0.44 


0.505 


0.041 


0.014 


0.016 


0.045 


0.44 


< 


x b 


< 


0.48 


U.OO ( 


U.U4Z 


U.Ulo 


U.Ulo 


U.U40 


0.48 


< 


x b 


< 


0.52 


0.677 


0.044 


0.016 


U.U11 


U.U5U 


0.52 


< 


x b 


< 


0.56 


0.796 


0.047 


0.017 


U.U3U 


U.U59 


0.56 


< 


x b 


< 


0.60 


0.991 


0.052 


0.018 


U.U56 


U.U79 


0.60 


< 


x b 


< 


0.64 


1.241 


0.058 


0.018 


U.U7U 


U.U92 


0.64 


< 


x b 


< 


0.68 


1.622 


0.068 


0.020 


U.U62 


U.U93 


0.68 


< 


x b 


< 


0.72 


2.092 


0.080 


0.028 


U.U44 


U.U96 


0.72 


< 


x b 


< 


0.76 


2.671 


0.094 


0.046 


U.U75 


U.128 


0.76 


< 


x b 


< 


0.80 


3.102 


0.104 


0.071 


U.14U 


U.189 


0.80 


< 


x b 


< 


0.84 


3.290 


0.111 


0.084 


U.2U1 


U.245 


0.84 


< 


x b 


< 


0.88 


2.953 


0.106 


0.065 


U.144 


U.19U 


0.88 


< 


x b 


< 


0.92 


1.897 


0.079 


0.094 


0.113 


U.167 


0.92 


< 


x b 


< 


0.96 


0.753 


0.042 


0.051 


U.2U5 


U.215 


0.96 


< 


x b 


< 


1.00 


0.090 


0.011 


0.004 


U.U61 


U.U63 



Table 6: The scaled _B-hadron energy distribution. 
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The statistical correlation matrix is shown in Table [7[ 
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-22.9 


-29.2 


-33.9 


-37.5 


-38.8 


-43.4 


-48.2 


-53.0 


-55.3 


-57.9 


-50.8 


-33.6 


-5.3 


37.5 


80.3 




21 


-15.3 


-22.3 


-27.3 


-30.8 


-33.8 


-35.7 


-40.3 


-44.2 


-47.4 


-50.1 


-57.7 


-64.0 


-62.5 


-51.3 


-19.8 


28.7 


77.8 


22 


-13.9 


-18.6 


-22.7 


-25.9 


-28.8 


-30.5 


-34.5 


-38.4 


-42.0 


-45.5 


-52.0 


-58.5 


-60.7 


-59.5 


-45.2 


-12.1 


39.2 85.6 


23 


-10.5 


-13.1 


-16.3 


-18.8 


-21.1 


-22.3 


-24.7 


-27.8 


-30.7 


-33.9 


-39.7 


-46.3 


-50.7 


-54.6 


-52.2 


-36.9 


0.8 52.1 85.8 


24 


-7.4 


-9.3 


-11.4 


-13.0 


-14.5 


-15.5 


-17.4 


-19.8 


-21.5 


-23.8 


-28.3 


-33.5 


-37.5 


-41.7 


-43.3 


-38.1 


-15.8 23.9 57.7 89.1 


25 


-5.1 


-6.4 


-7.7 


-8.1 


-9.0 


-10.3 


-11.8 


-13.7 


-14.7 


-16.0 


-19.2 


-22.9 


-25.8 


-28.7 


-30.9 


-30.3 


-19.0 6.2 32.9 68.8 91.2 



Table 7: The statistical correlation matrix (%). 
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7 Summary and Conclusions 



We have used the excellent tracking and vertexing capabilities of SLD to reconstruct the energies of B 
hadrons in e + e _ — > Z° events over the full kinematic range by applying a new kinematic technique to 
an inclusive sample of reconstructed 5-hadron decay vertices. The B selection efficiency of the method 
is 4.2% and the resolution on the B energy is about 9.6% for roughly 83% of the reconstructed decays. 
The energy resolution for low-energy B hadrons is significantly better than in previous measurements. 

We compared our measurement with several models of 6-quark fragmentation. The Bowler, Lund 
and Kartvelishivili et al. models, implemented within the JETSET string fragmentation scheme, de- 
scribe our data, as does the UCLA model. None of the Braaten et al, Collins-Spiller or Peterson et 
al. models implemented within JETSET, nor HERWIG, describes the data. 

The raw scaled 5-energy distribution was corrected for bin-to-bin migrations caused by the resolu- 
tion of the method, and for selection efficiency, to derive an estimate of the underlying true distribution 
for weakly-decaying B hadrons produced in Z° decays. Systematic uncertainties in the correction were 
evaluated and found to be significantly smaller than those of previous direct 5-energy measurements. 



The final corrected xb distribution D(xb) is shown in Figure jTTJ. This result is consistent with, and 



supersedes, our previous measurements |Tg, |17J . It is also consistent with a recent precise measure- 
ment 

It is conventional to evaluate the mean of this 5-energy distribution, < xb >■ For each of the 
seven parameter-dependent functions that provide a reasonable description of the data we evaluated 
< xb > from the distribution that corresponds to the optimised parameter (s); these are listed in 
Table |2| and Table [|. For the UCLA model, which contains no arbitrary parameters relating to b- 
quark fragmentation, we evaluated < Xb > from the corresponding unfolded distribution shown in 



Fig. [K§ this yields < xb > = 0.712. We took the average of the eight values of < xb > as our central 



value, and defined the model-dependent uncertainty to be the r.m.s. deviation. We obtained 

<x B > = 0.709 ± 0.003{stat.) ± 0.003{syst.) ±0.002 {model), (12) 

It can be seen that < xb > is relatively insensitive to the variety of allowed forms of the shape of the 
fragmentation function. 
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Figure 1: Distribution of the reconstructed P^-corrected vertex mass (points). The simulated distribu- 
tion is also shown (histogram) in which the flavor composition is indicated: b (open), c (cross hatched), 
and uds (dark shaded). 
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Figure 2: The relative deviation of the maximum missing mass from the true missing mass for simulated 
B hadron decays; the contributions from different B species are indicated separately: B° and B ± 
(open), B® (cross-hatched), and A& (dark shaded). 
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Figure 3: Distribution of the reconstructed MQ max for the selected vertices (points). The simulated 
distribution is also shown (histogram) in which the flavor composition is indicated: b (open), c (cross 
hatched), and uds (dark shaded). 
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Figure 4: Distribution of the reconstructed Ml max for the final selected sample (see text). The simu- 
lated distribution is also shown (histogram) in which the flavor composition is indicated: b (open), c 
(cross hatched), and uds (dark shaded). 
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Figure 5: The simulated efficiency for selecting B hadrons as a function of the true scaled 5-hadron 



energy, x% ue . 
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Figure 6: The fitted core and tail widths (see text) of the I?-energy resolution as a function of the true 
scaled I?-hadron energy. 
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Figure 7: Distribution of the reconstructed scaled .B-hadron energy (points) and the default Monte 
Carlo simulation (histogram). The solid histogram shows the simulated non-i? background. 
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Figure 8: The background-subtracted distribution of reconstructed 5-hadron energy (points). Also 
shown (histograms) in (a)-(f) are the predictions of the optimised models within JETSET (see text), 
(g), (h) show the predictions of HERWIG, and (i) of the UCLA model. Data points excluded from 
the fit are represented by open circles. 
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Figure 9: The background-subtracted distribution of reconstructed 5-hadron energy (points). Also 
shown (histograms) are the predictions of the optimised functional forms. Data points excluded from 
the fit are represented by open circles. 
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Figure 10: Distributions of the scaled weakly-decaying 5-hadron energy unfolded using different input 
models or functions (see text). 
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Figure 11: Final distribution of the weakly-decaying scaled B-hadron energy. In each bin the central 
value is the average of the eight distributions shown in Fig. the inner error bar represents the 
experimental error, and the outer error bar represents the sum in quadrature of the experimental and 
unfolding errors. 
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